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Abstract 

We review some of the properties of Advection-Dominated Accretion 
Flow (ADAF) models and show that they successfully describe many as- 
trophysical systems. Despite these successful applications some funda- 
mental problems still remain to be solved, the most important one being 
the physics of the transition between an ADAF and a geometrically thin 
Keplerian disc. 

PACS: 04.70.Bw; 97.60.Lf 
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1 Introduction 

'Advection Dominated Accretion Flows' (ADAFs) is a term describing infall of 
matter with angular momentum, in which radiation efficiency is very low. In 
their applications, ADAFs are supposed to describe inflows onto compact bodies, 
such as black holes or neutron stars; but very hot, optically thin flows are bad 
radiators in general so that, in principle, ADAFs are possible in other contexts. 
Of course in the vicinity of black holes or neutron stars, the virial (gravitational) 
temperature is Tvir « 5 x lO^^r"^ K, where r is the radial distance measured in 
the units of the Schwarzschild radius i?s = 2GM/c^, so that in optically thin 
plasmas, at such temperatures, both the coupling between ions and electrons 
and the efficiency of radiation processes are rather feeble. In such a situation, 
the thermal energy released in the flow by the viscosity, which drives accretion 
by removing angular momentum, is not going to be radiated away, but will 
be advected towards the compact body. If this compact body is a black hole, 
the heat will be lost forever, so that advection, in this case, acts like a 'global' 
cooling mechanism. In the case of infall onto a neutron star, the accreting 
matter lands on the star's surface and the (reprocessed) advected energy will be 
radiated away. There, advection may act only as a 'local' cooling mechanism. 
(One should keep in mind that, in general, advection may also be responsible 
for heating, depending on the sign of the temperature gradient - for example. 



1 



according to Nakamura et al. (1997), in some conditions, near the black hole, 
advection heats up electrons in a two-temperature ADAF). 

The role of advection in an accretion flow depends on the radiation efficiency, 
so that it depends on the microscopic state of matter and on the absence or 
presence of a magnetic field. If, for a given accretion rate, radiative cooling 
is not efficient, advection is necessarily dominant, assuming that a stationary 
solution is possible (see below). 



1.1 Simple ADAF solutions 

One can illustrate fundamental properties of ADAFs by a simple example. The 
advection 'cooling' (per unit surface) term in the energy equation can be written 
as 

Qadv = l^Csi (1) 

where M is the accretion rate, Cg the (adiabatic) speed of sound and ^ = 
— [(7 — l)~^(cilnT/(ilni?) — (dlnS/dlni?)] is the dimensionless advection fac- 
tor characterizing the entropy gradient (Abramowicz et al. 1995). 7 is the ratio 
of specific heats and we assume that the pressure is provided only by gas. 
Using the hydrostatic equilibrium equation 

I - - (2) 



where H is the disc semi-thickness and wk — \J GM/R = c/\/2r the Keplerian 
velocity, one can write the advection term as 

whereas the viscous heating term can be written as 

Qv. = ^n^ (4) 

where rh = (McKcs/47rGM) and Hcs is the electron-scattering opacity coefficient. 
Therefore, for very high temperatures, when (H/R) ~ 1 the advection term is 
comparable to the viscous term and cannot be neglected. 

From Eqs. (^ and (||) one can easily obtain an advection dominated solution 
by writing 

Qvis — Qadv (5) 

and using 

^m(aI])-V-i/^ (6) 



R J Kp 
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An ADAF is then described by the relation 

m = O.SSKesar-i/^i;. (7) 

In Eq. (^) we use the mass and angular momentum conservation equations and 
the viscosity prescription v = {2/?>)c^/ilK {^k = vk/R)- 

1.2 Maximum accretion rate 

Of course, the energy equation is 

Qvis — Qadv + Qrad (8) 

where Qrad is the radiative cooling (per unit surface), so that the existence 
of ADAF solutions depends on Qiad- There is no universal form for this term, 
which depends on the microscopic state of the matter, its content and the optical 
thickness. In the simplest case of an optically thin, one-temperature flow, cooled 
by non-relativistic free-free processes, 

Qrad - p'T'"^ - a-^m^-"" (9) 

so that the cooling term has a steeper dependence on m than both Qvis and Qadv 
One can see therefore that there exists a maximum value of m (Abramowicz et 
al. 1995) for which an ADAF solution is possible: m ^ a^r~^/^. We have 
assumed here that Eq. (^) applies for this value of m (see below). 

Abramowicz et al. (1995) found, in the non-relativistic free-free cooling case, 
that the maximum accretion rate for an ADAF solution is 

m,„ax = 1.7x 10Vr-i/2_ (iq) 

The value of mmax depends on the cooling in the flow and non-relativistic free- 
free cooling is not a realistic description of the emission in the vicinity (r 10^) 
of a black hole. The flow there will most probably form a two-temperature 
plasma (Narayan & Yi 1995). More realistic calculations in a 2T flow by Esin et 
al. (1997) give mmax ~ lOa^ with almost no dependence on radius. For larger 
radii mmax decreases with radius. 

1.3 Maximum a 

The existence of a maximum accretion rate for an ADAF is, as explained above, 
due to the form of the cooling law. For radii ^ 10^ optically thin cooling of 
gas pressure dominated plasma implies the existence of such a critical accretion 
rate. If, however, effects of optical thickness and of radiative pressure become 
important, the slope of Qrad(™) will change and, as shown by Chen et al. 
(1995), one may get families of solutions in which an ADAF is possible for all 
values of m because Qrad("i) has a maximum for Qrad < Qadv The existence 



3 



of such solutions depends on the value of a: rhmax exists only for a < amax, 
where the value of amax depends on the cooling mechanism and radius. As 
a simple example one can use a IT disc and consider that the simple free- 
free cooling formula is not valid for Tes ~ Ij where the effective optical depth 
Teff = y/ Kfi{Ks + Kes)'S,/2 and Kf£ is the free- free opacity. Using a Planck mean 
value for this optical depth, one obtains from rhmax < 'rh{Tf[ = 1) a condition 
for the viscosity parameter: 

a < acrit ~ r (11) 

Chen et al. (1995) obtained values of acrit between 0.2 and 0.4 for r = 5, 

depending on the model, and Bjornsson et al. (1996) get acrit > 1- If the 
viscous stress is assumed to be proportional to the gas, and not to the total, 
pressure, the value of acrit >> 1 (L0vas 1998) and is of no physical interest. 
It seems that for physically interesting configurations (a 1) there is always 
a mmax for ADAFs but one should not forget that the existence of such a 
maximum accretion rate is not a generic property of advective flows. 

1.4 Slim disc solutions 

For a < acrit there exists a separate branch of solutions which represents for m ^ 

1 the standard Shakura-Sunyaev discs. Higher accretion rates are represented 
by the so-called 'slim discs' (Abramowicz et al. 1988). It is sometimes said 
that slim discs represent solutions in which advection dominates because the 
optical depths are so high that photons are trapped in the flow. It is easy 
to see, however, that the optical depth of slim discs, for a rather large range 
of parameters, is rather low. In fact, it is the decrease of the optical depth 
with increasing accretion rate that is at the origin of the slim disc branch of 
solutions. Slim discs are only asymptotically advection dominated (see Fig 1. 
in Abramowicz et al. 1988). 

2 Global solutions 

The so called 'self-similar' solutions found by Narayan & Yi (1994) played a 
very important role in the development of ADAF astrophysics because they al- 
lowed Narayan and collaborators (sec Narayan, Mahadevan and Quataert 1998, 
for a recent review) to construct models which could be quickly compared with 
observations. These comparisons showed that ADAFs provide an excellent de- 
scription of such systems as the Galactic Center source Sgr A* (Narayan, Yi 
& Mahadevan 1995) and quiescent soft X-ray transients (SXTs) (Narayan, Mc- 
Clintock & Yi 1995). Of course, these models were rather crude and not self- 
consistent (for example, the 'advection parameter' /, which must be constant 
for a 'self-similar' solution varied with radius) but after more refined and con- 
sistent models had been calculated, they confirmed most of the prediction of 
the 'self-similar' ones. In fact, major revisions of the early models concerned 
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only the value of the transition radius between the ADAF and the outer geo- 
metrically thin disc, whose presence is required by observations of SXTs and 
the AGN NGC 4258 (see next section). 

The first 2T, global, optically thin ADAF solutions were obtained by Mat- 
sumoto, Kato & Fukue (1985). In Chen, Abramowicz & Lasota (1997) and 
Narayan, Kato & Fukue (1997) one finds the first detailed study of the prop- 
erties of global ADAF solutions. The Narayan et al. (1997) work deals with 
pure advection dominated solutions whereas Chen et al. (1997) consider more 
general solutions which allow for bremsstrahlung cooling. 

Figs. (|l|) and (||) show the characteristic properties of global ADAF solu- 
tions. They share some properties with slim disc solutions, such as the sub- 
Keplerian character of the flow for 'high' values of a and a super-Keplerian part 
of the flow for low values of this parameter. These features are related to the 
existence of a maximum pressure in the flow. Narayan et al. (1997a) argue 
that a pressure maximum is necessary for the existence of 'funnels' that appear 
in 'iron tori' (Rees et al. 1982). Since no flow models for iron tori exist it is 
difficult to tell the difference between these structures and ADAFs. 

Figs. (|l|) and (^) show also the influence of the maximum accretion rate for 
ADAFs: for low a's and/or high m's a flow which is advection dominated at 
small radii ceases to be an ADAF at larger radii. 

First studies of global ADAF solutions used the Paczyhski-Wiita (1980) 
approximation of the black effective potential, but such an approach is not 
satisfactory in the case of a rotating black hole. After the flrst ADAF solutions 
in a Kerr spacetime were found by Abramowicz et al. (1996), several other 
groups produced such ADAF models (Peitz & Appl 1997; Gammie & Popham 
1998; Popham & Gammie 1998). Jaroszyhski & Kurpiewski (1997) are the 
only authors who also calculated ADAF spectra in the framework of General 
Relativity. Fig (^) show Abramowicz et al. (1996) solutions for three values of 
the black hole angular momentum. One can see that these solutions are very 
similar to the "pseudo-relativistic" ones. In fact for a non-rotating black hole 
the Paczynski-Wiita (1980) ansatz is an excellent approximation. 

3 Applications 

Despite several rather serious uncertainties about the physics of ADAFs, their 
models are enjoying a growing fleld of application. The reason is simple. In many 
astrophysical systems powered by accretion (see below), the X-ray luminosity, 
which is supposed to probe accretion onto the central compact object, is very 
weak compared with the expectations based on independent estimates of the 
accretion rate and the assumption of high {'^ 0.1) radiative efficiency. The only 
model which predicts such properties of accreting systems is the ADAF model, 
of which a low radiative efficiency is a fundamental feature. 

Most models which reproduce observed spectra by ADAFs use the approach 
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pioneered by Ramesh Narayan and collaborators. It consists in fixing 'micro- 
scopic' parameters of the flow, such as the ratio of the magnetic to gas pres- 
sure, the fraction of viscous dissipation that goes directly to electrons, the a- 
parameter, the thermodynamical parameters, etc. Then the accretion rate is 
determined by adjusting to the observed X-ray flux. In a pure ADAF solution 
(no external disc) this determines the whole spectrum so that no additional 
freedom is allowed. The 'microscopic parameters' are considered to be universal 
so this procedure is, in fact, an one-parameter 'fit' (assuming of course that 
the black hole mass is given by independent observations). In the case of a 
two-component flow (ADAF -|- accretion disc) a second 'parameter' has to be 
flxed: the transition radius between the two flows. In principle this should not 
be a free parameter, but should be given by the physics of this transition. This 
is however a weak (the weakest according to the present author) point of the 
ADAF 'paradigm': the mechanism of transition between the two flows is un- 
known. There exist the 'usTial suspects', such as evaporation, but a consistent 
physical model is yet to be found. A principle according to which the flow will 
be an ADAF when, at a given radius, an ADAF and an accretion disc solutions 
are possible has been used (see e.g. Menou, Narayan & Lasota 1998) but the 
case of NGC 4258 seems to contradict this principle (see Sect. 3.3). 

3.1 Sgr A* 

Rees (1982) was the first to suggest that accretion onto the Galactic Center 
black hole might be advection-dominated. Detailed models by Narayan et al. 
(1995,1998) and Manmoto, Mineshige fc Kusunose (1997) provide rather im- 
pressive fits to observations from radio to X-ray frequencies. In this case the 
model is a pure ADAF; no outer accretion disc is present. 

3.2 Soft X-ray Transients and 'Low state X-ray binaries' 

Soft X-ray Transients (SXTs) are a natural field for application of ADAF models 

since an accretion disc model is unable to fit the observations of these systems in 
quiescence (Lasota 1996). In the first ADAF model of quiescent SXTs proposed 
by Narayan, McClintock & Yi (1996) an outer stationary accretion disc was 
responsible for the observed emission in optical and UV frequencies. It was 
shown however, by Lasota, Narayan & Yi (1996), that it is impossible to find 
an ADAF -|- accretion disc configuration which would be consistent both with 
observations and with the requirements of the disc instability model, which is 
supposed to describe SXT outbursts. A new version of the model was proposed 
by Narayan, Barret & Yi (1997) in which the contribution of the outer disc was 
negligible. This was achieved by increasing the transition radius and increasing 
the magnetic to gas pressure ratio so that the optical/UV emission is due to 
synchrotron radiation from the ADAF. 
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The validity of this model has been independently confirmed by Hameury 
et al. (1997). They reproduced the multi-wavelength properties of the rise to 
outburst of GRO J1655-40, including the observed 6 day delay between the 
rise in optical and X-rays, by using the disc instability code of Hameury et al. 
(1998). As initial conditions they used parameters given by an ADAF fit to the 
quiescent spectrum. They also showed that it is very difficult (if not impossible) 
to reproduce these observations with different initial conditions. 

Narayan (1996) proposed to interpret various spectral and luminosity states 
observed in X-ray binaries using ADAF + accretion disc models with varying 
accretion rate and transition radius. This idea has been applied to the black 
hole SXT X-ray Nova Muscae 1991 (Esin et al. 1997) and to Cyg X-1 (Esin 
et al. 1998; the first to apply an ADAF model to this system was Ichimaru 
1984). The existence of a two-component flow with varying accretion rate and 
transition radius in the XNova Muscae 1991 was confirmed by observations of 
the X-ray reflected component in the Ginga spectra of this object, but these 
observations suggest (Zycki, Done & Smith 1998) that the Esin et al. (1997) 
model requires some modifications. 

3.3 NGC 4258 

This LINER is a very important testing ground for the ADAF models. The 
black hole mass in this system (3.6 x 10 ''Mq) is very well determined due to the 
presence of narrow, water maser lines (Miyoshi et al. 1997). The observed X-ray 
luminosity is ~ 10~^LEddington and the bolometric luminosity is no more than an 
order of magnitude larger. Lasota et al. (1996b) proposed, therefore, that the 
inner accretion in this system proceeds through an ADAF. New observations 
in infrared (Chary & Becklin 1997) and a new upper limit on the radio flux 
(Herrnstein et al. 1998) constrain the transition radius to be rtr ~ 30 (Gammie, 
Blanford & Narayan 1998). In this model the accretion rate would be m = 
9 X 10""^, much higher than the value of to ^ 10^^ proposed by Neufeld & 
Maloney (1995) in their model of the masing disc, but in agreement with values 
obtained by Maoz & McKee (1997) and Kumar (1997) for the same disc. It 
seems that the value proposed by Neufeld & Maloney (1995) is excluded by the 
IR and X-ray data. In NGC 4258, an ADAF solution is possible for r > rtr so 
that the validity of the principle according to which an ADAF is preferred over 
an accretion disc whenever the two solutions are possible seems questionable. 

3.4 LINERs and weak AGNs 

Lasota et al. (1996b) suggested that also other LINERS and 'weak' AGNs 

could contain ADAFs in their inner accretion regions. A recent analysis of the 
variability of such systems (Ptak et al. 1998) brings new arguments in favour 
of this hypothesis. 
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3.5 Nuclei of giant elliptical galaxies 

Fabian & Rees (1995) proposed that the weh known problems with understand- 
ing the properties of nuclei of giant elliptical galaxies, which emit much fewer 
X-rays than one would expect from independent estimates of the accretion rate, 
could be solved if the accretion flow formed an ADAF. Reynolds et al. (1996) 
and di Matteo & Fabian (1997) applied ADAF models to the nuclei of M87 
and M60 respectively. The latter authors point out that observations of the 
predicted spectral turnover at ~ 3 x 10^^ Hz puts constrains on the values of 
the 'microscopic' parameter of the ADAF. 

3.6 BL Lac's 

In BL Lac's the emission is dominated by the jet but several properties of these 
systems points out to the possible presence of an ADAF (Madejski 1996; Celotti, 
Fabian & Rees 1997). 

3.7 Remnant and MACHO black holes in the ISM 

Fujita et al. (1998) noticed that if black holes accreting from the interstellar 
medium formed an ADAF, the characteristic ADAF spectral form would help 
in their detection. 

4 Conclusion 

ADAF models find successful applications in many domains of accretion astro- 
physics. They are the only models which describe both the dynamics of the 
accretion fiow and its emission properties. In any case, all models which require 
the presence of optically thin, very hot plasmas (such as 'coronal' models) must 
take into account advective heat transport in order to be self-consistent (see 
Section 1.2). ADAFs are the only solutions that satisfy this requirement. 

Of course, as mentioned above, there are still serious problems to be solved 
before one can conclude that ADAFs, in their present form, are the models of ac- 
cretion flows around black holes at low accretion rates. The main problem to be 
solved is that of the transition between an accretion disc and an ADAF. Further- 
more, in order to reproduce data, ADAF should form two-temperature plasmas 
with ions at virial temperature and much cooler electrons. This requires weak 
coupling between ions and electrons and implies that viscosity heats mainly the 
ions. Some recent studies suggest that it is difficult to achieve this together with 
equipartition between magnetic and gas pressures (Quataert & Gruzinov 1998). 
Such calculations however, concern extremely complicated physical processes so 
that it is not clear that results of mmierical simulations apply to real systems. 
One should therefore apply a more pragmatic attitude, that the apparently suc- 
cessful applications of ADAFs to many astrophysical systems suggests that no 
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mechanisms coupling ion and electron exist in real astrophysical configurations 
(Fabian & Rees 1995). 

The validity of ADAF models will be decided by observations. 
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Figure 1: Disk solutions for a = 0.1 and m = 10~^ - solid lines and a = 0.1 and 
m = 10~^ - dashed lines. The mass of the black hole is lOM©. M. is the Mach 
number, f7 the angular frequency of the flow and ^ the advcction parameter. 
The heavy lines are the global solution and the thin lines are the corresponding 
self-similar solution. Note that M., O, and ^ depend on m very weakly. 
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Figure 2: Solutions for a = 0.01 and m = 10^^ - solid lines and a — 0.001 
and TO — 10^^ - dashed lines. The black hole mass is the same as in Fig. 
(0). Heavy lines correspond to global solutions and thin lines correspond to 
self-similar solutions. Note the super-Keplerian angular momentum and the 
maximum pressure near the transonic region. Note also that for a = 0.001 and 
rh = 10~^, the local cooling becomes important for large radii. 
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Figure 3: The radial structure of the Mach number {A4), the pressure {p in 

cgs unit), the angular momentum {C in unit of M), and the sound speed (cg 
in unit of the speed of light). The mass of the black hole is IOMq, a = 0.1 
and M /Me = 10~^. The solid, dotted, and dashed lines represent the cases 
of a/M = 0, 0.5, and 0.99 respectively. The heavy dots represent solutions 
obtained with the pseudo-Newtonian potential. These solutions are excellent 
approximation to the solutions representing Schwarzschild black hole flows (in 
the case a = 0) . The corresponding Kcplcrian angular momenta of test particles 
around Kerr black holes are also shown for comparison (the thin lines). 
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